Boson features in STM spectra of cuprate superconductors: Weak-coupling phenomenology 
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We derive the shape of the high-energy features due to a weakly coupled boson in cuprate superconductors, 
as seen experimentally in BiiSr-zG 'aiCu20s+x (BSCCO) by Lee etal [Nature 442, 546 (2006)]. A simplified 
model is used of d-wave Bogoliubov quasiparticles coupled to Einstein oscillators with a momentum indepen- 
dent electron-boson coupling and an approximate analytic fitting form is derived, which allows us (a) to extract 
the boson mode's frequency, and b) to estimate the electron-boson coupling strength. 
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Scanning tunneling microscopy (STM), applied to the su- 
perconducting cuprate I^S^CaiC^Os+a; (BSCCO 2212) 
yl, found a feature in the density of states (DOS) at an en- 
ergy well above the energy scale of the so-called coherence 
peak energy (Fig. [TJ, and attributed it to an electron-boson 
coupling. In conventional (s-wave) superconductors (e.g. Hg, 
Pb, Al), such features due to electron-phonon coupling were 
known in tunneling spectra from superconductor-insulator- 
normal metal junctions Q 0]. The phonon frequencies in- 
ferred from the tunneling feature agreed with the phonon den- 
sity of states inferred from neutron scattering; furthermore, 
the phonon-mediated superconducting T c and gap were cor- 
rectly predicted [2] from the tunneling using the Eliashberg 
formalism fl. In the case of cuprates, the mechanism for 
superconductivity is not established, and there are divergent 
opinions whether the mode observed by Lee et al contributes 
to the pairing |d 0-H • 

In BSCCO, the pairing strength is highly inhomogeneous at 
the nanoscale 111 llll2ll . as inferred from the spatial fluctuations 
of the energy _E co h of the "coherence peak" in STM spec- 
tra (Figure [T). Lee et al discovered that the boson feature's 
energy Et, os "floats" with the same inhomogeneity as £? co h, 
namely i?bos = £-coh + f&lo with a (spatially uniform) boson 
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FIG. 1: [COLOR ONLINE] A typical measured STM spectrum in 
BSCCO (proportional to DOS n{uj)) as a function of energy E; the 
data was provided by Jacob Alldredge. Energies of the "coherence 
peak" -Ecoh and boson feature iJbos are indicated. Boxes show en- 
ergy windows used to fit the analytic form (see Fig. [3] below.) 



frequency flo- To infer i?bos, they identified it as the inflection 
point in n(u>) before the feature. In this paper, we improve on 
this recipe by deriving an analytic formula for the boson fea- 
ture, starting from the simplest phenomenological model of 
a cuprate and using basic RPA calculations. Our focus here 
is the energy dependence rather than the spatial modulations 
H3. [7^7 of this feature. One prior calculation jH] addressed 
the same question of extracting Hq from from the shape of 
the DOS of BSCCO, using more elaborate (Eliashberg) cal- 
culation, but in an entirely numerical framework, making the 
interpretation harder and the method impractical for extensive 
fitting. 

We first ask just what point in the feature is to be identified 
as Et, os : our recipe implies a value for Ml in agreement with 
the analysis in Ref. . Secondly, we ask how can one extract 
the electron-boson coupling strength; our results indicate it is 
indeed small enough that our weak-coupling approximation 
is justified, and furthermore this coupling alone is unlikely to 
explain the magnitude of the superconducting gap. 

Weak-coupling Model — 

We begin by setting up the simplest possible model, tak- 
ing the electron-boson coupling as a small perturbation to 
an already superconducting fermion dispersion of the stan- 
dard mean-field form (as in Ref. 1 131) . and then setting up the 
DOS calculation within the RPA approximation. Our anal- 
ysis is agnostic as to the boson's nature, which is sometimes 
argued to be magnetic (9D, hut usually considered to be an 
oxygen vibration, on account of the O 18 isotope effect Jl]]. 

Our bare fermion Hamiltonian has the usual mean-field 
form 



H = ^ e(k)<4. a c Ka + A(k)c k , CT c_ k ^ 



h. 



(1) 



k.cr 



where e(k) is the normal-state band dispersion, for which 
(in all numerical calculations in this paper) we adopt Nor- 
man's six-parameter tight-binding fit to ARPES data on 
BSCCO 11 1 CHI . The quasiparticle dispersion is then E(k) = 
\J e(k) 2 + A(k) 2 , where we will assume e?-wave pairing with 



A(k) ^ ^ [cos(M - cos(fc„)] . 



(2) 



We (plausibly) approximate the bosonic mode as a disper- 
sionless (Einstein) oscillator at frequency £lo, and assume an 
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electron-phonon coupling 



rte-ph — 



7= 5(q)4+q, CT C k, CT (&-q + &q) 



(3) 



k,q,(7 



where 6 q and 6 q are the bosonic creation and annihilation op- 
erators, and N is the number of lattice sites. For simplic- 
ity we work through the case <7(q) = g; after completing 
that, we will revisit the more general case with a momentum- 
dependent p(q). 

Our object, the DOS, is defined as the trace of the electron 
term in the Green's function: 



n(u) = --Tr k ImGn(k,w), 



(4) 



where Trk = fl 2 / BZ ePk/(2-7r) 2 , and the integral is over the 
Brillouin zone. In the 2x2 Nambu formalism, the bare 
Green's function is given by 



£°(^)-=(- (k) 



e(k) A(k) 



to + e(k) 



(5) 



We shall henceforth use Pauli matrices r it and adopt the gauge 
in which A is real: thusG = [M+e(k)r 3 +A(k)r 1 ]/[w 2 - 
E(k) 2 )]. The boson propagator has the form 



i 



= £>(fi). (6) 



Self-Energy and density of states due to the boson — 
The boson feature enters the DOS via the dressed Green's 
function, in the RPA approximation, 



Gtk,^)- 1 = G^k,^)- 1 -E(k,w). 



(7) 



Because Oq an d D(Q.) were momentum-independent, so is 
the electronic self energy, reducing (at lowest order in g) to 
£(k, lu) = where 

SH = .9 2 J ^^(«)Tr q {r 3 G (k-q; W -r!)T3} (8) 

After a contour evaluation of the integral Eq. (H) reduces 



to 



EM 



-Tr k 



(u + O )I + e(k)T3 - A(k)r x 



(w + fi ) 2 - £(k) 2 



B(k) 



AM 
B(k) 



2j) 



[ W -£(k)] 2 -fi 2 



(9) 



The off-diagonal (tj) terms in Eq. (O vanish, E12 = E21 = 
0, since Ak has <i-wave symmetry (reverses sign under 90° 
rotations). 

We write n(u) — no(ui) + 5n(uj), where no(u>) is tbe basic 
DOS in the absence of the phonon coupling [derived from ©I 
and has the well-known "coherence peaks" centered at energy 
values i-Ecoh close to ±Aq; Sn(uj) contains contributions of 
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FIG. 2: [COLOR ONLINE] Numerically computed DOS as a func- 
tion of energy. The standard six-parameter dispersion for BSCCO 
was used [10]; we chose the pairing amplitude to be Ao = 0.2|ti | « 
29.6 meV and the boson energy to be Mlo = 0.25jti | ~ 37 meV, 
where t\ is the nearest-neighbor hopping. Three different choices of 
damping are shown: from sharpest to flattest, g = 0.005, 0.01 and 
0.02 (in units of ti). Lower inset shows real and imaginary parts 
of the electron self-energy function. This has singularities at i?bos 
rounded by the damping. 



order g 2 , in particular the boson feature. Writing the Taylor 
expansion of (0, we extract the terms in G\\ linear in E and 
thus 



8n(uj) 



-— ImTri 

7T 



- e (k)] 2 E 11 H + |A(k)| 2 S 22 H - 

T O T7T/1 \919 



= — Im{/iMEiiM + J2ME22M} • 

7T 



(10) 



This is the first version of our result, suitable for numerical 
fits 01511 . but requiring integrations over the zone at each inter- 
ation [for the key formulas (O and ( TTOb . Note that in numeri- 
cal calculations, we replace lu —> oj + ig in ©, where ig rep- 
resents the physical quasiparticle damping (from all sources 
except our boson mode), a parameter found essential for fit- 
ting the "coherence peaks" in the DOS [16fl. (It is easy to 
replace this energy-independent damping by r)(cu), as used in 
Ref. iflql .) Fig.|2]shows a representative numerical calculation 
of the self-energy function (inset) and the resulting DOS. We 
see a dip-hump shape, in agreement with experiment; -Ebos 
falls between the dip and the hump similar to the assumption 
of Ref. El. 

Asymptotic Form near E\ ms = -E co h + f2n — 
We now extend our results to an approximate analytic for- 
mula, for the boson feature's shape, by treating not only the 
electron-phonon coupling g, but also the damping 77 as a small 
parameter: in the limit 77 — ^ the feature is a singularity cen- 
tered at E hos = Eco h + Q . 

First recollect the origin of the familiar "coherence peak" in 
the basic DOS no(cu): it is a van Hove singularity due to the 
saddle points at k = (k S) n) and equivalent momenta where 
the Fermi surface crosses the zone boundary. The pertinent 
pole in G is ^[lj — E^k)^ 1 !; there is no contribution from 
t 3 due to the factor e(k) which vanishes on the Fermi sur- 
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face. It is well known that Trk[w — £(k)] at a saddle gives a 
logarithmic singularity, so we find a singular part 



-Re/(w - E coh + i?7coh) (1 1) 



with 



f(z) = Hm*z/4K x K y ). 



(12) 



Here E(k) » £ coh + (k x - fc s ) 2 /2 mi - - 7r/a) 2 /2r 



near the saddle, to* 



i, x ni.y, and i^, X y are cut-offs, 



v /TOim ?/ 

representing the range of (k x , fc y ) within which this expansion 
is valid. For our parameters, 1/m* = 94.07meV a 2 , and we 
take K x = 0.5a -1 and K y = 0.06a _1 for later numerical 
calculations. 

The self-energy E(-E) has a singularity due to the same sad- 
dle point, with the pole of form {\) 2 g 2 l[u - E(k) ~ f2o]~\ 
coming from the second big term in Clearly, integrating 
over k gives the same logarithmic divergence, with its argu- 
ment shifted by flo. Thus, 

■ 2 2 

S(w+i?7) = regular terms +— — m* f(u>—E\, os +irf)I (13) 

with f(z) from (TTZb . This behavior is confirmed by the inset 
of Fig. 

The I dependence in ( fT3l signifies that En ss £22 at the 
singularity. Thus ( fTOb simplifies to 



5n sin s( w ) = -Im[£ u (w)J(£ boB + Mj)], 



(14) 



with 7(cj) = [7i (w) + 1 2 (w)/2 [see Eq. [TO) 

Thus our key asymptotic result is that <5n(u;) has a logarith- 
mic singularity at i?bos, rounded by the finite damping r\. The 
result is a linear combination of a rounded step and a cut-off 
log divergence, with the exact shape (and the location of Et, os 
within it) depending on the phase angle in = \I{uj)\e t 'l' 1 , 
which depends on the band structure [cf. Eq. 1141 1. 

At u) = 115 meV, the rough dependence on damping is 
I(uj + irj) = 1.5 x 1(T 5 (?/ - 15) + 0.7 x lCT 3 i. Thus, the 
shape of Sn(ui) is a (comparable) combination of a rounded 
upwards step from ReSn and a rounded logarithmic hump 
from ImSn leading to location of the boson mode frequency 
u) before the hump (as seen in numerics cf. Fig. |2|. 

We can attach physical interpretations to the real and imag- 
inary parts of £22(1^). The imaginary part represents an in- 
elastic event in which a real phonon excitation is created; the 
real part represents the quasiparticle being dressed by virtual 
phonons. 

Since the predicted feature includes a "step up", we are in 
agreement with the recipe of Lee et al which placed Et, os at 
the inflection point before the hump of the boson feature, mo- 
tivated by previous work on molecular vibrational features in 
electron tunneling |18, 19]. Refs.0 and0 located Ebos even 
lower, at the minimum of the dip in the dip-hump feature. As 
mentioned before, we also place Et, os before the hump but 
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TABLE I: Fit parameters for the "coherence peak" using Eq. l |15t 
(left column) and for the boson feature using Eq. l !16b (right col- 
umn). 



more specifically in between the hump and its preceding in- 
flection point. 

We can attempt to compare our self-energy functions with 
those of Ref. [5] [(Figure 3(c)], computed numerically from 
Eliashberg theory. ReSii(w) is proportional to their ImZ(w) 
which indeed resembles a (positive) log divergence, while 
ImEii oc 1 — Z(lu) shows a rounded up step. 

Fitting Scheme for the Experimental Boson Feature — 

In this section, we translate our asymptotic forms to a sim- 
plified fitting scheme for our weak-coupling model and, by 
applying it to the experimental spectrum in Fig. Q] extract 
the Ehos and also obtaining the electron-phonon coupling g 
from the boson feature's amplitude. 02111 We consider the 
experimental signal to be in arbitrary units so we write it 
h(uj) = f3 ca in(ui), where the coefficient f3 ca \ includes un- 
known factors such as the STM tip set-point. As the dis- 
persion e(k) is already known from ARPES [10], the "coher- 
ence peak" is sufficiently constrained that we can calibrate 
/3 ca i from it. We read off E co ^ = 40.8 meV from the peak 
position in Fig. Q] From this, using E co h = E(k sa ddic), we 
infer A = 44.23meV. 

The saddle point of the quasiparticle dispersion at k sa( jdie 
contributes a logarithmic singularity to the DOS at the "co- 
herence peak": 



n (u) 



(15) 



where Uq ng (w) is given by (TTTb . and we adopt the simplest us- 
able form n rcg (w) = a co hW + b co h for the regular part, which 
is due mainly to no(w). 

Table [Ogives the results of the calibration fit to the data in 
Fig. Q] using energies in (30 meV, 50 meV). As Fig. [3] (left 
panel) shows, the fitting is good in this window. This fit gives 
a quasiparticle broadening 77 co h sa 10 meV (assumed to be 
constant over the Brillouin zone and the energy window 30 — 
50meV), uncomfortably large in that f] co h/E co h m 1/4. We 
do not know why this exceeds the result 7?(E co h w 40meV) « 
1 meV. fitted by Ref. 16 assuming a broadening r/(w) oc ui. 

Now we turn to the fit of the boson feature, using an energy 
window (80 meV, 140 meV) which contains the hump in Fig. 
Q] to the fitting form implicit in Eqs. ( fTOb [for I(u>)], ( fT3l . and 
(O: 

™M = <o g sM + ^ sins M- (16) 

Here we take the simplest usable form for the regular part 



4 



3 




30 



a) 



£coh=40.8 
oj (meV) 



50 



b) 



£ bos =96.41 
oj (meV) 



135 



FIG. 3: [COLOR ONLINE] Fit of the experimental DOS n(w) to 
fitting forms, using the windows of energies marked in Fig.Q] (a) Fit 
of "coherence peak" to Eq. j 1 5b . using energies 30-50 meV. (b) Fit 
of boson feature to Eq. dl6t . using energies 80-140 meV. 



n^Jui) — abos^ + febos, representing no (w) plus regular con- 



tributions from Also from (fl4l we see 



2ig 2 a 2 7 



xlm 



(2^) 4 

[I(u + if] hos ) ■ f(w - E hos + iVbos)] 



(17) 



The fitted parameters are given in Table [0 the fit (Fig. O is 
fairly good in its energy window. 

We note that, based on the data from which Fig.[T]is drawn, 
Ref. id identified the bosonic mode energy as 52 ± 8 meV, 
using the inflection point before the hump, so our result of 
~ 56 meV (fitting just one typical spectrum) is in agreement 
with them. The quasiparticle damping was ?7bos ~ llmeV. 
Thus rj/E(k) s» ri/Eb os s» 0.11 <C 1 in the E^ fit window, 
verifying the criterion for the Bogoliubov quasiparticles to be 
well-defined. 

A dimensionless measure of coupling strength is the the 
ratio of the logarithmic factors f(z) in the boson feature 
[Eq. ( [ToT il and coherence peak [Eq. (fl5T ll: 



Ai, 



_ 2 3 2 |/(£, 



OS 



ir)bc 



(2tt) 2 



0.057, 



(18) 



using the numerical value \I(E\ )os + i?7bos)| = 8.7 
meV -2 , validating our weak-coupling assumption. 



x 10" 



Is^)! 2 — 9 2 [^( cos 1x + cosq y )], where we set g to the fitted 
g value from TableU We compute the gap renormalization us- 
ing the obvious generalization to account for q dependence of 
electron-boson coupling in (the off-diagonal components of) 
Eq. [8] We find that for all energies, the gap renormalization 
is less than 5 meV, which is small enough compared to Ao to 
justify our weak-coupling assumption, but not so small to cat- 
egorically rule out some contribution by the boson to pairing. 

For the boson feature, the overall structure of the calcula- 
tion carries through but the self-energy £ becomes momen- 
tum dependent. We find the same sort of DOS feature, in 
which "g 2 " is now interpreted as a certain weighted aver- 
age of |g(q)| 2 over the Brillouin zone - a lumped parameter 
in the spirit of the "a 2 F(uj)" combination from the strong- 
coupling formalism J2]. The singularity in the self en- 
ergy still come from the saddle point in the dispersion of the 
(i-wave BCS quasiparticles, leading to the same qualitative 
shape (smoothed step + logarithm) for the boson feature. 

Conclusion and Discussion — 

We have shown how a weak-coupling point of view can be 
used to analyze the high-energy features in the STM data of 
BSCCO. The ideal analytic shape of the feature is a linear 
combination of a (rounded) logarithmic-kink and a (rounded) 
step edge [cf. Eq. (fl4l l. Our proposed fitting scheme al- 
lowed us to extract (1) the boson's frequency 51o (2) an aver- 
age electron-boson coupling g, and an estimate of the damp- 
ing of the d-wave Bogoliubov quasiparticles. Our estimate 
Hq m 56 meV is in agreement with previous estimates from 
STM data, which were not fully in agreement with ARPES 
data. Our simplified simple functional form for the boson fea- 
ture [Eq. ( fToT ll facilitates the vast number of numerical fits re- 
quired by the extreme spatial inhomogeneity of STM spectra 
in BSCCO iGjIGKiei]. However, our theory did not address 
the spatial Fourier spectrum of the boson feature Ill [l3i[l4ll . 
which might distinguish the true functional form of <?(q) and 
thus illuminate the nature of the bosonic mode. 

In an s-wave superconductor with a dispersionless boson 
such as ours, the boson feature would cancel completely in 
our weak-coupling treatment. This may have motivated sug- 
gestions ||9l|2fj(] that the observed phonon mode was an artifact 
localized in a surface barrier. However, we found the d-wave 
quasiparticle dispersion is non-uniform enough to make S(o>) 
nonconstant in cj thus producing an DOS feature. As noted 
after, Eq. (0 so long as we have a momentum-independent 
electron-phonon coupling to an Einstein oscillator, no off- 
diagonal self-energy contribution is generated, and hence the 
d-wave gap is not renormalized. 

Our approach was agnostic as to the pairing mechanism. 



Momentum- dependent boson coupling and gap renormalizationlf the fitted g respects the weak-coupling assumption, as we 

it can be inferred that the boson producing the STM 



What if the electron-boson coupling g(q) in ([3]l is not 
constant but depends on the electron momentum transfer q? 
Firstly, it gives renormalizations of A(k) due to £12 which is 
no longer zero (see Eqs.(|4|) and (01). To obtain an upper bound 
for the gap renormalization, we try the form for q depen- 
dence which leads to the maximal renormalization, namely 



found, 

feature is not contributing to the pairing; if the weak-coupling 
assumption is violated, we can only conclude that the boson 
perhaps plays a role in the mechanism. To resolve that ques- 
tion, one must see if a proper strong-coupling Eliashberg cal- 
culation predicts a pairing amplitude Ao comparable to the 
observed value. 
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